Geometry of fully coordinated, two-dimensional percolation 
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We study the geometry of the critical clusters in fully co- 
ordinated percolation on the square lattice. By Monte Carlo 
simulations (static exponents) and normal mode analysis (dy- 
namic exponents), we find that this problem is in the same 
universality class with ordinary percolation statically but not 
so dynamically. We show that there are large differences in 
the number and distribution of the interior sites between the 
two problems which may account for the different dynamic 
nature. 
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I. INTRODUCTION 

The geometrical phase transition known as percolation 
(see, for a review, Stauffer and Aharony ^]) is appreci- 
ated by many to be an elegant and simply defined yet 
fully featured example of a second order phase transi- 
tion. A number of variations of the original percolation 
problem were proposed as better models of some physical 
phenomena in the past. This includes the backbone perco- 
lation for studying electrical conduction through random 
media, polychromatic percolation for multi-component 
composites, and four- coordinated bond percolation for 
hydrogen-bonded water molecules. In particular, Blum- 
berg et al II and Gonzalez and Reynolds § studied a 
random bond, site-correlated percolation problem they 
call four-coordinated percolation on the square lattice. 
They conclude that this problem belongs to the same uni- 
versality class as the ordinary random percolation with 
the same set of (static) exponents. 

In this paper, we revisit a problem in this realm, 
though not exactly the same one. We define fully coordi- 
nated percolation as the site percolation problem where 
only the occupied sites all of whose neighboring sites are 
also occupied can transmit connectivity. Since the ran- 
dom element is the site, this problem is slightly different 
from the bond problem referred to above. Thus, after 
generating a random site configuration with the indepen- 
dent site occupation probability p, we only select those 
occupied sites with all 4 neighbors also occupied on the 
square lattice and study the clusters formed by nearest 
neighbor connections among those sites. It should be 
noted that this problem is distinct from the so-called 
bootstrap percolation (see, e.g., |J) where sites of less 
connectivity are iteratively removed. In our problem, 
no iterative procedures are involved; rather, sites of less 



than full connectivity are marked first and then all of 
them removed at one time. 

This problem arose in the context of studying the vi- 
brational properties of fractal structures tethered at their 
boundaries In that problem, scaling was observed 

in the normal mode spectrum whose origin may lie in the 
ratio of 2 length scales, one of which is the size of highly 
connected regions of a cluster. In this context, we have 
embarked on revisiting the characteristics of randomly 
generated, but highly connected geometrical structures. 

In the next section, we summarize the Monte Carlo and 
finite size scaling analyses of the static critical properties 
of fully coordinated percolation. In Section 3, we discuss 
the normal modes of the transition probability matrix for 
tracer diffusion on the structure using the methods of 
Arnoldi and Saad (see, e.g., ||). Then in Section 4, we 
describe the classification of the cluster sites into external 
boundary, internal boundary, and interior ones and using 
these to show the major distinctions between the critical 
clusters of ordinary and fully coordinated percolation. 
We summarize the results in the final section. 



II. STATIC CRITICAL BEHAVIOR 

To determine the static critical behavior of fully coor- 
dinated percolation we first performed Monte Carlo sim- 
ulations on a square lattice in two dimensions. Each site 
is occupied with probability p independently and subse- 
quent fully coordinated sites are marked and their con- 
nectivity searched. Lattice sizes of I? where L = 256, 
512, 1024, and 2048 were constructed. For each lattice 
size we further made a thousand realizations wherein a 
different random number seed was used on every run. 

The unnormalized susceptibilities, i.e., S(L) = '^ JS s' 2 h s 
where n s is the number of clusters of size s, are calcu- 
lated on each run and are then summed at the end of 
the thousand realizations. The average susceptibilities 
X are calculated by dividing the sum by the number of 
realizations and the lattice size. The prime on the sum- 
mation indicates the fact that the contribution of the 
largest cluster to x near and above what we perceived 
to be the critical probability p c has been subtracted as 
usual ||. 

In Fig. 1 we plot the average susceptibilities against 
the probability p for the corresponding lattice sizes. The 
data correspond to the values of L = 256, 512, 1024, and 
2048 from the lowest to highest. We can see that the 
effects due to the finite sizes of the lattices are exhibited 
clearly. In particular, there are well-defined peaks which 
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scale with lattice sizes as 

X(Pmax,L) ~ Utl" (1) 

where the known exact value of 7/1/ for the ordinary per- 
colation is || ~ 1.7917. To demonstrate the precision 
of our calculations, we plot xiPmax, L) against the corre- 
sponding lattice sizes in the inset of Fig. 1. Notice that 
the data follow an excellent power law, leading to a least 
squares fit of xiPmax, L) ~ L 17911 . The value of 7 found 
is identical with the ordinary percolation value to within 
about 0.03%. This result confirms previous work jH|] 
stating that fully coordinated percolation and ordinary 
percolation belong to the same static universality class. 

The critical behavior of susceptibility is known to scale 
as § 

X(p,oo) ~ \p - pep 7 (2) 

where for ordinary percolation 7o = f§ ~ 2.3889. Notice 
however that in Fig. 1 the peaks are very near p = 1.0. 
This would provide data to the right of the peaks in only 
a small probability interval. In our simulations, we would 
therefore use x only to the left of the peaks. 

Since the scaling relation in Eq. is expected only 
for infinite lattices, we use only the data taken from L = 
2048 to test it. Since there are two unknowns in Eq. (||), 
we first choose a particular p c and make a fit to see what 
value of -f exp is obtained. If we choose p c — 0.886 we get 
lexp = 2.4004. The correlation coefficient, \R\, for this 
fit is 0.99999. The discrepancy between ^y exp and -f Q is 
around 0.481%. Choosing p c = 0.8858 we obtain j exp = 
2.3864. The discrepancy this time is around 0.10% and 
\R\ = 1. So we have an exact fit for this value of p c and 
the jexp found is very close to the 7 . Choosing p c — 
0.885 we obtain j exp = 2.3302 with an \R\ = 0.99999. 
The discrepancy for this value of p c is 2.46%. Fits done 
with p c between 0.885 and 0.886 gave \R\ = 1; however, 
the j exp found when p c = 0.8858 gave the closest value to 
7 . This allows us to conclude that 7 for fully coordinated 
percolation is the same as that for ordinary percolation 
while also giving an estimate for the value of p c close to 
0.8858. (We will state the experimental uncertainty for 
p c after all our analyses are presented.) 

From the fit done to examine the scaling in Eq. (|IJ) we 
could further conclude that v for fully coordinated per- 
colation should be the same with that for ordinary per- 
colation. This again confirms the statement that fully 
coordinated percolation is in the same static universality 
class as ordinary percolation. Another universal constant 
often used to characterize ordinary percolation is the am- 
plitude ratio C+/C- of susceptibility \ (whose value is 
about 200 in d = 2 0). In fully coordinated percolation, 
this quantity is unfortunately difficult to calculate accu- 
rately because the critical region for p > p c is very small 
(see below). When we constrain the exponent 7 to be 
close to 7 and use the p c estimated in this work, how- 
ever, we find that C+ /C_ is of i9(10 2 ), which is consistent 
with the above observation as well. 



The contribution of the largest cluster to the suscepti- 
bility is not significant when p < p c . However when p ~ 
Pc a significant number of sites will belong to the largest 
cluster and when p > p c the largest cluster is dominant in 
the whole lattice. The average susceptibility contribution 
due to this largest cluster is xi = 2 s Li/(^ 2 ^)i where 
the summation is over ./V = 1000 realizations and s max 
is the size of the largest cluster. The fractal dimension, 
df, can be obtained from s max by 

SmaxiPc) ~ L d f (3) 

where Smax(Pc) is the mean size of the largest cluster at 
p c - Xi should therefore scale as 

Xi ~ L 2d f - \ (4) 

For ordinary percolation on a two dimensional lattice 
(see, e.g., ||), it is known that df — 91/48 and y = 
2df — 2 ~ 1.7917. For fully coordinated percolation, the 
scaling in Eq. (Q) have two unknowns, p c and df. Simi- 
lar to what we have done when examining the scaling in 
Eq. (^|) , we choose trial values for p c and then perform a 
least squares fit to obtain the corresponding y = 2df — 2. 
By looking for the range of trial p c that maximizes the 
regression coefficient \R\, we arrive at an estimate of p c 
to be close to 0.8845 where \R\ = 1 and y = 1.7855. 
The variation of \R\ is about 2 parts in 10 5 if p c is var- 
ied by 0.0002, always with less than 1% deviation from 
the ordinary percolation value of y. From these results 
we conclude that df for fully coordinated percolation is 
the same as that for ordinary percolation as well as the 
estimate of about 0.8845 for p c . 

In addition to the above, we have also performed the 
scaling analysis of the quantity %i(p, L), as both p and 
L are varied, in the form of 

Xn = L 2d f - 2 g(\p - PcV L) (5) 

where g(\p — p c \ u L) is a scaling function. Using the 
exactly known ordinary percolation values of the expo- 
nents df and v (as they have been shown to be the 
same for fully coordinated percolation above), we ob- 
tained the maximum data collapsing in the range of 
0.884 < pc < 0.885. 

Independent of the above analyses based on the fully 
coordinated clusters obtained by Monte Carlo simula- 
tions of fixed-sized square grids, we have also performed 
Monte Carlo simulations by growing fully coordinated 
clusters starting from a seed site using a variant of the 
breadth-first search algorithm This latter approach 
has an advantage that there is no obvious finite-size ef- 
fects and that statistics taken while a cluster is still 
growing represents a partial sum automatically. That 
is, we start growing such clusters 10,000 times at each of 
p = 0.880, 0.884, 0.885, and 0.890, and keep track of how 
many of them are still growing at predetermined inter- 
vals of size (2™ where n = 1, 2, ... , 15 in our case). This 
number, say, N s represents the partial sum 
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TV./iVi = - y s'n s , (6) 

y s'>s 

Since the normalized number of size-s clusters, n s 
scales as s~ T f(es a ) where t = 187/91 and a — 36/91, 
we expect that N s scales as 

N s s T - 2 ~ f(es°) (7) 

near p c and for large s. In particular, at p c , this quan- 
tity should be constant independent of (large) s. The 
numerical results are shown in Fig. 2, where the data 
correspond, from highest to lowest, to p = 0.890, 0.885, 
0.884, and 0.880. The horizontal dashed line drawn to 
guide the eye makes it clear that data for p = 0.885 best 
approximates a horizontal line as s — > oo, suggesting that 
a good estimate of p c would be 0.885. 

We now consider all the above results together. The re- 
sults from scaling in Eq. (g) indicate a range 0.885—0.886, 
and those from scaling in Eq. (Q) indicate 0.8844-0.8847, 
while those from scaling in Eq. (|5|) hints p c to be in the 
0.884 — 0.885 interval. Another result that could also be 
used are the values of p for the peaks in Fig. 1, which 
vary from 0.8841 to 0.8844 (with the peak for the largest 
grid L — 2048 occurring at p pea k — 0.8844). Combining 
all these results, our final estimate is p c = 0.885 ± 0.001. 

III. DYNAMIC CRITICAL BEHAVIOR 

By dynamic critical behavior here we simply mean the 
asymptotic long-time behavior of diffusion taking place 
on an incipient infinite cluster of fully coordinated per- 
colation, or its equivalent scalar elastic behavior. This 
represents the simplest kind of dynamics associated with 
these complicated geometrical objects and is mainly re- 
flected in the two dynamic critical exponents called d s 
(spectral dimension) and d w (walk dimension). 

It is well known that the return-to-the starting point 
probability of the random walk, P{t), in the long-time 
limit obeys the power law, 

P(t)~t- d °/ 2 , (8) 

where d s is the spectral dimension of the walk. In a frac- 
tal medium, d s is less than the space dimension d, be- 
cause the progressive displacement of the random walker 
further from the starting point is hampered by its en- 
counter with the irregularities of the medium at all scales. 
Thus, d s is expected to be greater for environments that 
provide higher connectivity at large length scales, inde- 
pendently of the fractal dimension itself which is mainly 
the measure of the overall size scales or how many sites 
are connected, not how well those sites are connected to 
each other. 

For media with long-range loops, df (fractal dimen- 
sion), d s , and d w are not independent but are expected 
to obey the well-known Alexander and Orbach scaling 
law @ 



d s = 2d f /d w . (9) 

For this reason, we only calculate d s here though both 
of d s and d w can be conveniently calculated by numeri- 
cally studying the transition probability matrix W which 
represent the random walk on a specific fractal medium. 
Our calculation in this work is only one aspect of such 
an analysis: finite size scaling of the dominant non-trivial 
eigenvalue which describes the longest finite time scale of 
the Brownian process. This approach has already been 
described in detail elsewhere [|| , and thus we merely state 
the main feature and then immediately report our spe- 
cific numerical results. 

The matrix W is constructed from the elements 
being equal to a hopping probability per step (equal to 
■| here) for available nearest neighbor sites i and j. For 
each neighbor site which is not present, a probability of 
\ is added to the probability for not taking a step for 
one time period - this is called the blind ant rule. Many 
large matrices W are obtained by Monte Carlo simula- 
tion (by growing a fully coordinated percolation cluster 
from a seed site and stopping the growth when a prede- 
termined desired size is reached) and their largest eigen- 
values are numerically obtained by the so-called Arnoldi- 
Saad method. The dominant non-trivial eigenvalue Ai 
is the largest eigenvalue just below the stationary eigen- 
value 1 and it is known to satisfy the following finite size 
scaling law: 

|lnA 1 |wl-Ai-5'- 2/<i3 . (10) 

Shown in Fig. 3 are our results from such an analy- 
sis. We have generated at least 1000 independent real- 
izations of the underlying fully coordinated percolation 
clusters for sizes S = 1250, 2500, 5000, 10000, and 20000 
at each of the three nominal probabilities p — 0.883, 
0.885, and 0.887, and numerically obtained Ai for each 
cluster. The main part of Fig. 3 shows the data from 
p = 0.885, the value shown to be closest to p c in this 
work. The figure shows an excellent power law fit (re- 
gression coefficient of -0.99993) to Eq. (|To|) with the ex- 
ponent 2/d s = 1.486 ± 0.01. This power translates 
to d s = 1.346 ± 0.011 which is close to but definitely 
larger than the corresponding ordinary percolation value 
of di°^ — 1.30 ± 0.02 estimated by many independent 
calculations. (See, e.g., Q.) For comparison, a loopless 
variant of percolation has exactly the same static expo- 
nents as ordinary percolation but has d s w 1.22 in two 
dimensions, about twice as much deviation from ordinary 
percolation in the opposite direction as the present fully 
coordinated percolation problem. jll) 

In the inset for Fig. 3, we show a normalized Ai by 
plotting (l—\i)S 2 / dl ° } where the circles are for p — 0.887, 
squares for p = 0.885, diamonds for p = 0.883 and the 
solid line is a horizontal line (for ordinary percolation) 
to guide the eye. In all cases, the standard errors of 
the mean for each set of data are substantially smaller 
than the size of the symbols used in the figure. The 
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distribution of Ai in each case appears to be Gaussian 
with the standard deviations scaling in the same way as 
the means. 

^From these results, we conclude that, though the nu- 
merical differences are small, there is a high likelihood 
that the fully coordinated percolation clusters are sig- 
nificantly different from the ordinary percolation coun- 
terparts even at long length scales. In the next section, 
we show that this analysis is vindicated by exposing one 
dramatic difference in the cluster morphology which will 
not be obvious to an uncritical observer. 



IV. CLUSTER GEOMETRY 

In this section we examine the geometry of fully coordi- 
nated percolation clusters more closely. First, we present 
Fig. 4 which show in grey scale the sites of (a) fully coor- 
dinated percolation cluster and (b) ordinary percolation 
cluster at respective p c . The overall visual impression is 
that they are very similarly shaped even down to the de- 
tails of the boundaries and internal holes. Their shapes 
are also essentially independent of the underlying lattice 
anisotropy (as Fig. 4(a) has actually been rotated by 45 
degrees with respect to the coordinate axes of the square 
lattice). However, the number and distributions of the 
especially dark points are evidently quite distinct in (a) 
and (b). They cluster more and are much more abundant 
in (a) than in (b) . These sites are actually the interior or 
fully coordinated sites in the internal part of the cluster. 
The remaining sites (shaded grey) are either the external 
hull sites or internal boundary sites. 

In Fig. 5 quantitative examination is made on the dif- 
ferent classes of sites of the two kinds of clusters. In 
the main part of the figure, the average numbers of two 
kinds of sites are shown, interior (diamonds for fully co- 
ordinated percolation, crosses for ordinary percolation), 
and external (squares for fully coordinated percolation, 
plusses for ordinary percolation). It is clear that the in- 
terior sites are more than 3 times abundant in fully co- 
ordinated percolation as is visually suggested by Fig. 4. 
Though this is primarily a local effect due to the full co- 
ordination rule, they do have a multiplicative effect at 
long-range connectivity and thus may well be the source 
of the small difference in the value of d s . 

Of course just the fact that there are more than 3 
times as many interior sites (and correspondingly, much 
fewer hull sites) in fully coordinated percolation must 
have quantitative consequences (even if not qualitative) 
for any process on the cluster which depend on degrees of 
connectivity rather than just on the number of connected 
sites. An example of the effect of the different numbers 
of the hull sites may be in oxidation or catalysis of a ma- 
terial through the external embedding phase or even a 
irregularly shaped breakwater in the form of the exter- 
nal boundary of a percolation cluster. The external sites 
are those which are sometimes called hull sites and they 



are known to scale with an exact exponent in ordinary 
percolation as 

N hull ~ S d */ d f (11) 

where — 7/4 and thus the exponent is 84/91 = 

0. 923... for ordinary percolation. Since this is less than 

1, these sites comprise less and less fraction of the clus- 
ter as s — > oo and the remaining sites (i.e., interior and 
internal boundary sites) eventually dominate the whole 
cluster. This is already evident from the greater slopes 
close to 1 for the interior sites in Fig. 5. The linear 
regression fits for the hull sites in Fig. 5 indicate the 
slopes of about 0.913 for ordinary percolation and 0.922 
for fully coordinated percolation with essentially perfect 
fits, again reinforcing the conclusion that they show the 
same static critical behavior. 



V. SUMMARY AND CONCLUSION 

In summary, we have studied both static and dynamic 
critical behaviors associated with a model of the highly 
connected regions of a disordered cluster. The model is 
a site variant of the four-coordinated percolation HU on 
the square lattice we call fully coordinated percolation. 
While the bond version was studied for static critical be- 
havior, neither bond nor site version was previously stud- 
ied for the dynamic behavior to the best of our knowl- 
edge. We have used various methods such as Monte Carlo 
simulations, finite-size scaling and Arnoldi-Saad approx- 
imate diagonalization of large random matrices for this 
purpose. 

Though all indications are that the static behavior of 
this model is exactly the same as the ordinary percola- 
tion (as previous work suggested), the dynamic behavior 
shows a small but significant difference in the values of 
the universal critical exponents. We have looked for the 
cause of this difference and found a three-fold increase in 
the number and significantly enhanced clustering of the 
interior sites (i.e., those not on the exterior or internal 
boundaries) and the associated decrease in the number of 
boundary sites. Thus, although the deviations from or- 
dinary percolation in terms of the values of the dynamic 
critical exponents are not large, there will be rather sig- 
nificant differences in any processes that depend sensi- 
tively on those numbers. Possible examples of such pro- 
cesses include the oxidation of a material through the 
external embedding phase and the vibrational normal 
modes with boundary conditions such as clamping or 
tethering of the external boundaries (through the con- 
trast in elastic constants of embedding and embedded 
materials, for example). 
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FIG. 3. The scaling of the largest nontrivial eigenvalue Ai 
of the Markov chain analysis is shown by plotting 1 — Ai 
against the cluster size s on a log-log scale. The inset shows a 
normalized A (see text) for p = 0.887, 0.885, and 0.883 (from 
top). The solid line is horizontal, corresponding to the slope 
for ordinary percolation. 
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FIG. 4. (a) A typical fully coordinated percolation cluster 
at its p c , and (b) the corresponding critical ordinary perco- 
lation cluster. The dark points are the interior sites of full 
coordination and the grey shaded sites are either on the ex- 
ternal hull or internal boundaries. 

FIG. 5. The numbers of external sites (squares and plusses) 
and those of interior sites (diamonds and crosses) are scaled 
against the cluster size s. The solid (fully coordinated perco- 
lation) and dashed (ordinary percolation) lines are the linear 
least squares fits. In the inset, we show the number of ex- 

d <o) /d <0) 

tcrnal sites divided by s h f , showing that the population 
of the external hull sites of fully coordinated and ordinary 
percolation scale in the same way. 



FIG. 1. Susceptibilities \ for fully coordinated percolation 
on finite square lattices of size L 2 are shown against proba- 
bility p, where crosses, diamonds, squares, and circles are for 
L = 256, 512, 1024, and 2048, respectively. The inset shows 
the finite-size scaling of the peaks of these curves by plotting 
the maximum x vs - L on a log- log scale. 



FIG. 2. The scaled partial sum N 3 is plotted against s a for 
p = 0.890 (cross), 0.885 (circle), 0.884 (triangle), and 0.880 
(plus). The dashed line is a horizontal line to guide the eye, 
showing that the data for p = 0.885 best fits the horizontal 
slope for large cluster sizes. 
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Fig. 1 (Cuansing, Kim & Nakanishi) 




Fig. 3 (Cuansing, Kim & Nakanishi) 



